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Abstract 

We present a general class of spatio-temporal stochastic processes describing the causal 
evolution of a positive-valued field in space and time. The field construction is based on 
independently scattered random measures of Levy type whose weighted amplitudes are inte¬ 
grated within a causality cone. General n-point correlations are derived in closed form. As 
a special case of the general framework, we consider a causal multiscaling process in space 
and time in more detail. The latter is derived from, and completely specified by, power-law 
two-point correlations, and gives rise to scaling behaviour of both purely temporal and spa¬ 
tial higher-order correlations. We further establish the connection to classical multifractality 
and prove the multifractal nature of the coarse-grained field amplitude. 
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Nonspecialist summary 


Throwing many dice many times at different points in space is an ex¬ 
ample of an uncorrelated random proeess, because the number of points 
on a particular die is independent of those on other dice around it, and 
because the dice do not have any memory. Such uncorrelated random 
processes may seem unsuitable for describing correlated phenomena in 
nature. Levy-based modelling, however, accomplishes exactly that. It 
does so by making use of overlapping sums of dice as follows: Suppose 
we have three dice labeled A, B and C whose outcomes are uncorrelated. 
The two variables X = A + B and Y = B + C will nevertheless be correlated 
since X and Y share the outcome of die B. This simple example can be 
generalised to construct correlated processes in spacetime. Let, for exam¬ 
ple, the energy e at a point (x, t) be determined by the sum of outcomes 
of all dice occurring within its ambit set, a kind of causality cone similar 
to Einstein’s familiar light cone (see Figure 1). As shown in Figure 2, 
two energies at different points will then be correlated if their respective 
ambit sets overlap, because they will share a common ancestry to some 
extent. The freedom to choose both the kind of randomness and the form 
and size of the ambit set permits this approach to mimic different types of 
correlated behaviour. This article deals with the particular class of phe¬ 
nomena called multifractal and multiscaling, which includes turbulence, 
data traffic flows, cloud distributions, rain fields and tumour profiles, to 
name but a few. All of these show a power-law-like behaviour in their 
correlations which are easily incorporated into the Levy-based modelling 
scheme by taking products of random variables rather than their sums. 
We show how to construct both the kind of randomness (the Levy ba¬ 
sis) and suitable ambit sets for this class, taking as input the measured 
correlations, and calculate analytically correlations for overlaps of various 
kinds (see Figures 2 and 3). While this paper concerns itself chiefly with 
multifractals, the Levy scheme as such is much more general and can be 
appled to many other correlated random processes. 


1 Introduction 

Multifractality ^ has in the last decade become one of a number of well-established approaches 
to the analysis of time series and spatial patterns, whether nonlinear, random, deterministic, or 
chaotic. It serves, for example, to characterize the intermittent fluctuations observed in fully 
developed turbulent flows HE] and in data traffic flows of communication networks H ■ Spatial 
patterns of cloud distributions and rain helds reveal multifractal properties, as do super¬ 
rough tumor profiles [7|. While still controversial, multiscaling has also been applied to financial 
time series of exchange rates and stock indices lainiiini. Many other examples may be found 
in the literature. 

Multifractality should not, however, be seen as a mere tool for analysis and characterization: 
it has also found its way into theoretical modeling. Maybe the simplest construction of a 
multifractal field is achieved with random multiplicative cascade processes H which introduce 
a hierarchy of scales and multiplicatively redistribute a flux density from large to small scales. 

Various generalizations of such purely spatial and discrete cascade processes towards continu¬ 
ous cascade processes in time and/or space, formulated in terms of integrals over an uncorrelated 
noise field, have been undertaken recently. A purely temporal and causal generalization to a 
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continuous cascade process is, for example, discussed by Schmitt m, who introduces a log¬ 
normal field, itself defined as an integral of a weighted and uncorrelated noise field over an 
associated time-dependent interval. By judicious choice of integration interval and weight func¬ 
tion, the resulting process is stationary and exhibits approximate scaling behaviour of two-point 
correlations. 

Muzy and Bacry m discuss a similar approach, constructing a purely temporal multifractal 
measure with the help of a limiting process. Since, however, the field amplitude depends on 
times later than the observation time t, the model does not obey causality. 

The proper description and modeling of spatio-temporal multifractal physical processes 
clearly calls for a model generalization that is causal, explicitly depends on space and time 
and does so in a continuous framework. A first step in this direction was achieved in m, where 
a continuous and causal spatio-temporal process was constructed in analogy to a discrete cas¬ 
cade process. Analytical forms for two- and three-point correlations for the case of a stable noise 
field were successfully compared to the corresponding experimental statistics in fully developed 
turbulent shear flow. 

The aim of the work presented here is to provide a general framework for the construction 
of spatio-temporal processes that permits a unified description of the above-mentioned models 
[niiiadsi while transcending them all. The basic notion in this framework is that of inde¬ 
pendently scattered random measures of Levy type. The appealing mathematics behind these 
measures, as described in [Tlj (with emphasis on spatio-temporal modeling), provide a charac¬ 
terisation of arbitrary n-point correlations independent of the choice of a concrete realisation of 
the model. This opens up the possibility of designing spatio-temporal processes almost to order, 
i.e. satisfying prescribed correlations. 

As an application, we present the construction of a multiscaling and causal spatio-temporal 
process that is based on and derived from scaling two-point correlations. In contrast to m and 
m, where the specification of the probability density of the noise-field must be included from 
the very beginning, we can construct the process without fixing the marginal distribution of 
the field-amplitude. This opens up the possibility of tailoring the marginal distribution of the 
process to the phenomenology of a given application. In particular, the special case of a stable 
law coincides with m- 

While we will concentrate on multifractal examples in most of this paper, it should be noted 
that this framework is not restricted to multiscaling (defined as scaling of correlation functions) 
or multifractal processes (defined as scaling of the coarse-grained process). 

The paper is structured as follows. In Section [2 we discuss the general framework for 
spatio-temporal modeling and derive an explicit expression for n-point correlation functions for 
the general set-up. Based on this result, we turn to the application in the context of causal and 
multiscaling spatio-temporal processes in Section 13 where we show in detail the multiscaling 
properties of temporal and spatial n-point correlations of arbitrary order and establish a relation 
between spatio-temporal multiscaling and spatio-temporal multifractality. Section 0] concludes 
the paper with a summary and a brief outlook. 

2 General model approach 

The aim of this Section is to define the general framework and to provide useful mathematics for 
the construction of a class of causal spatio-temporal processes that are based on the integration 
of an independently scattered random measure of Levy type. The integral constituting a given 
observable extends over a finite domain in space-time, called the ambit set S. This approach 
includes the special case of a continuous cascade process in space and/or time as an example. In 
particular, we recover the temporal cascade processes discussed in m and m, as well as the 
spatio-temporal cascade process derived in m In this paper, we restrict ourselves to causal 
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processes in 1 + 1 dimensions only, referring the reader to m and na for the general case of 
(n+l)-dimensional processes and its various applications and properties. 

The basic notion is that of an independently scattered random measure (i.s.r.m) on continous 
space-time, RxR. Loosely speaking, the measure associates a random number with any subset 
of R xR. Whenever two subsets are disjoint, the associated measures are independent, and 
the measure of a disjoint union of sets almost certainly equals the sum of the measures of the 
individual sets. For a mathematically more rigorous definition of i.s.r.m.’s and their theory of 
integration, see Refs. [nmniiizi. 

Independently scattered random measures provide a natural basis for describing uncorrelated 
noise processes in space and time. A special class of i.s.r.m.’s is that of homogeneous Levy bases, 
where the distribution of the measure of each set is infinitely divisible and does not depend on 
the location of the subset. In this case, it is easy to handle integrals with respect to the Levy 
basis using the well-known Levy-Khintchine and Levy-Ito representations for Levy processes. 
Here, we state the result and point to m for greater detail and rigour. 

Let Z be a homogeneous Levy basis on RxR, i.e. Z(S) is infinitely divisible for any S C 
R X R. Then we have the fundamental relation 


exp 


h{a)Z{da] 


= exp 


K[/i(a)]da > , 


( 1 ) 


where (• • •) denotes the expectation, h is any integrable deterministic function, and K denotes 
the cumulant function of Z{da), defined by 


In (exp {^Z{da)}) = K[^] da. 


( 2 ) 


The usefulness of o is obvious: it permits explicit calculation of the correlation function of the 
integrated and /i-weighted noise field Z{da) once the cumulant function K of /i is known. 

2.1 General model ansatz 

Based on relation o, we construct a spatio-temporal process that is causal and continuous^ by 
defining the observable field e{x,t) as 


e{x, t) = exp 


I S{x,t) 


h{x, t] x', t')Z{dx'xdt') 


(3) 


This is clearly a multiplicative process of independent factors exp{L(x, t] x', t') Z{dx' xdt')} made 
up of a specifiable weight function h and a homogeneous Levy basis Z over RxR. Contributions 
to field amplitude e{x,t) lie within the influence domain S{x,t), called the associated ambit set. 
To guarantee causality, we demand that S be nonzero only for times preceding the observation 
time t, i.e. S{x,t) C Rx[—oo,t] (see Figure 1). 

Ansatz reduces to the model of Ref. m when focusing on one purely temporal dimension, 
setting S{t) = [t -|- 1 — A, t], A > 1, /i(t; t') = {t — and defining Z to be Brownian motion. 

Similarly, a non-causal and again purely temporal version of the general model with a conical 
ambit set leads^ to the scale-dependent measures used in m- 

As shown in the Appendix and discussed in Section mi our approach also includes the case 
of a multifractal measure that is constructed without a limit-argument. Moreover, it allows for 
multifractality in space and time simultaneously. This multifractal case (with the additional 
assumption of a stable Levy basis) corresponds to the log-stable process described in where 

^In this context, continuity refers to the definition of observable t{x,t) for a continuous range of points {x,t). 

^This connection can be established by replacing the spatial coordinate x with a scale label and omitting the 
causality condition. 
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the ambit set is constructed from an analogy to a cascade process. In Section HOI we derive the 
same result from an alternative approach. 

Some other applications of © are discussed in m. 

The generality of the model is based on the possibility of choosing the constituents of 
the process e(x, t) independently. The available degrees of freedom are the weight function h, 
an arbitrary infinitely divisible distribution for the Levy basis Z (including Brownian motion, 
stable processes, self-decomposable processes etc.) and the shape of the ambit set S. As all of 
these quantities can be chosen to fit the purpose and application in mind, our ansatz permits 
sensitive and flexible modeling of the correlation structure of e{x,t). Despite its generality, the 
model is tractable enough to yield explicit expressions for arbitrary n-point correlations in closed 
form. 


2.2 n-point correlations 

The dehnition of the process e(x, t) allows for an explicit calculation of arbitrary spatio-temporal 
n-point correlations, defined as 


Cni^l ) ) ■ • • ) Xm tn) — (®1 j ) ' ■ • • ' tn}) j 


which give a complete characterisation of the correlation structure of e(x,t). 
Using the definition ® and the fundamental relation o we rewrite 

Cn (2^1, , . . . , Xn ,tji) — 


= (exply^ [ h{xi,ti‘,x',t')Z{dx'xdt')\\ 

\ ]/ 

Z(dx'xdto|^ 


= ( exp 


= exp 


K 


/RxR 


^ ^ ^S{xi,ti)h{Xi, ti, X , t ) 


'^i=l 


dx'dt' > , 


where we made use of the index-function 


\A{x,t) = 


1 when (x, t) G A 
0 otherwise 


(4) 


(5) 


( 6 ) 


for sets A C RxR. The last step in © follows from the fundamental equation ©• 

To illustrate ©, we consider in more detail the cases n = 2 and n = 3 with the abbreviation 
Si = S{xi,ti)- For n = 2, it follows that 


(e(xi,ti)e(x 2 ,t 2 )) = exp < / K [/i(xi, ti; x, t)] dx dt 

[JSl\S2 J 


X exp ■ 


'S2\Sl 


K [h{x 2 ,t 2 ; X, t)] dx dt 


X exp < / K[h{xi,ti; x,t) + h{x 2 ,t 2 ', x,t)] <ix dt> . 

USinS2 J 


(7) 


As illustrated in Figure 2.c, the first and second factor are contributions from the non-overlapping 
parts of the ambit sets, while the third stems from the overlap of S{xi,ti) and S{x 2 ,t 2 ) (the 
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shaded area). The latter factor describes the correlation of the field amplitude e at different 
spatio-temporal locations; for locations where the overlap 5 (xi,ti) fl S{x2,t2) vanishes, we get 
uncorrelated field amplitudes {e{xi,ti)e{x2,t2)) — (e(xi, ti))(e(x2,^2)) = 0 . Thus, the extension 
and shape of S{x,t) characterises the range of correlations, while the overlap of two ambits, the 
weight-function h and the cumulant function K influences the correlation strength. 

In third order we get a similar result. The combinatorics of the overlap of ambit sets 
for the observation points (xi,ti), (x2,t2) and (xsjta) yields seven disjoint domains as fol¬ 
lows: the three domains S{xi,ti)\[S{x2,t2) U S{x3,t3)], S{x2,t2)\[S{xi,ti) U 5(x3,t3)] and 
S{x3,ts)\[S{xi,ti) U S{x2,t2)] give uncorrelated contributions associated solely with one field 
amplitude; for instance, S{xi,ti)\[S{x2,t2)U S{x3,t3)] is the contribution to e{xi,ti) that is in¬ 
dependent of e(x2, ^2) and e(x3, ts). A second set of three domains [S{xi,ti)r\S{x2,t2)]\S{x3,t3), 
[ 5 (xi,ti) r\ S{x3,t3)]\S{x2,t2) and [S{x2,t2) r\ S{x3,t3)]\S{xi,ti) constitute the contributions 
to the correlation of two field amplitudes but without that of the third field amplitude. Finally, 
S{xi,ti) n S{x2,t2) n S{x3,t3) is the overlap of all three ambit sets that describes the common 
correlation of all three field amplitudes. 

Using the simplified notation = K[h{xi-^ , ; x, t)+h{xi2 ,ti2-,x,t)-\ - \-h{xi ., U .; x, t)], 

the result in third order hence reads 

{e{xi,ti)e{x2,t2)e{x3,t3)) = 


= exp < / Ki dx dt > exp 

[JSl\iS2US3) ) 


/ K2 dx dt 

/S2\(SlUS3) , 


exp < / K3 dx dt 

[JS 3 \iSlUS 2 ) 


X exp 


< / Ki^2dxdt > exp < / 


X exp I f Ki^2,3dxdt| . (8) 

USinS2nS'3 J 

It is clear from the above examples that the correlation structure corresponds directly to an 
intuitive geometrical picture in which the design and the overlap of the ambit sets S determine 
the correlation structure. 

Conversely, one can use some given correlation structure Cn as the starting point for designing 
a suitable shape of the ambit set and weight-function h to fit these requirements, opening up a 
wide range of applications. As outlined in the next section, multiscaling appears as a specific 
example, while Ref. m provides further insight into the kind of processes that can be modeled 
and explores the potential of the additive counterpart defined as lne(x,t)). 

3 Multiscaling model specifications 

In this Section, the concept of multiscaling and multifractality is examined in the context of the 
general model approach presented above. An explicit expression for the ambit set S is derived 
from scaling two-point correlations, and fusion rules m expressing n-point correlations solely 
in terms of scaling relations are formulated. Finally, the link to standard multifractality is 
established. 

3.1 General remarks and assnmptions 

In order to keep the mathematics as transparent as possible, we will use some simplifying as¬ 
sumptions about the structure of the process e(x, t). Our goal is the construction of a stationary 
and translationally invariant process with scaling two-point correlations. For the simplest way 
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to achieve stationarity and translational invariance, we assume h = 1 and take the form of the 
ambit set S{x,t) to be independent of the location {x,t), so that 

S{x,t) = {x,t) + So , (9) 

where the shape of Sq is independent of {x,t). (Note that /i = 1 is not a prerequisite for 
stationarity and translational invariance. It would be sufficient to require h{x', t'] x, t) = h{x', t'), 
but we can do without this additional degree of freedom for the special case of scaling relations 
for two-point correlations.) 

Figure 1 illustrates the various features of S{x, t), which we now discuss. At the origin (0,0), 
it is specified mathematically by 

So = {{x, t) € RxR : —T < f < 0, —g{t + T) < x < g{t + T)} . (10) 


This definition contains a finite decorrelation time T, ensuring that no correlation survives for 
temporal separations At larger than T, e.g. (e(x,t)e(x,t -|- At)) — (e(x, t))(e(x, t -|- At)) = 0 for 
all At > T. 

Spatially, the ambit Sq is limited by a function g{t), whose monotonicity ensures that the 
spatial extension of the causality domain increases monotonically for past times. The non¬ 
constancy of g implies a time-dependent spatial decorrelation length /(At), since, when two 
observations are separated by a space-time distance (Ax, At) (as illustrated in Figure 2.c), the 
two-point correlation (e(x,t)e(x -|- Ax,t -|- At)) — (e(x,t))(e(x -|- Ax,t + At)) vanishes for all 
Ax > /(At) = g{At) + g{0). The spatial decorrelation length /(At) decreases monotonically 
with At, and its maximum /(O) = 2g{0) = L defines the decorrelation length L. This is a 
physically desirable property. 

Finally, we impose a locality condition g{T) = 0, i.e. the ambit set Sq is attached to (x, t) in 
an unequivocal way. 

The procedure followed in the next section starts from the assumption that spatial and tem¬ 
poral two-point correlations scale, and constructs the model according to this requirement. The 
basic relation we use in the translationally invariant and stationary case under the assumptions 
@ and /i = 1 is 


(e(x, t)e{x + Ax, t + At)) 



= exp I [ K[l]dxdtl exp I [ K[l]dxdtl exp | [ (K[2] — 2K[l])dxdf 

US'! J U52 J USinS2 


= (e)2exp|y(Ax,At)(K[2]-2K[l])|, (11) 

where we have used 0 with h = 1 and the abbreviations Si = S{x,t), S 2 = S{x + Ax, t -|- At) 
and 

V (Ax, At) = Vol(5'(x, t) n 5(x -|- Ax, t -|- At)) (12) 

for the Euclidean volume of the overlap of the ambit sets. Due to translational invariance and 
stationarity, we have (e(x,t)) = (e(x -|- Ax,t -|- At)) = (e). 

Assuming that K[2] > 2K[1] (note that, by the strict convexity of log-Laplace transforms, 
we always have K[2] — 2K[1] > 0), Eq. (II Ijl can be solved for V, 


V{Ax,At) 


K[2] - 2K[1] 


In 


(e(x, t)e(x -|- Ax, t -|- At)) 

w 


(13) 
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This relation establishes a simple geometrical way to design a model with prescribed two-point 
correlations: one has only to choose ambit sets S{x,t) in a way that the volume of the overlap 

fulfils (lllfl) . This will be done in the next Section for the case of scaling two-point correlations 

(see m for more examples other than scaling relations). 

3.2 Construction of the ambit set via scaling two-point correlations 

Implementing the general framework © together with the above assumptions and procedure, we 
start out by demanding power-law scaling for the lowest-order spatial and temporal correlations, 

{e{x,t)e{x + Ax,t)) = CxiAx)~'^^‘^\ Ax € [heal, L^cai] C [0,L], (14) 

(e(x, t)e{x, t + At)) = Ci(At)"^(^\ At G [tscauTscai] C [0, T], (15) 

with Cx and ct constants. Note that the scaling exponents r(2) appearing in (ITll) and m are 
taken to be identical; differing spatial and temporal scaling exponents, as used previously e.g. 
in m, are easily accommodated within our model, but do not satisfy the simpler relations m 
given below. 

Following the recipe sketched in m, we get, using stationarity, for the temporal two-point 
correlation m the expression (see Figures 1 and 2.a) 

F(0, At) = / 2g{t)dt = / 2g{t)dt -|- / 2g{t)dt 

J lAt J lAt '^T^~T'scal 

_ Inct — ln((e)^) r(2)lnAt 

K[2] - 2K[1] ~ K[2] - 2K[1] 

for At G [tscai,Tscai], and after differentiation of both sides with respect to At, we 
expression 

^ 2(K[2] -2K[1]) t’ * ^ [hcal,Tscal] 

for the function g{t) bounding the ambit set S{x, t) within the temporal scaling regime [tscai,Tscai]- 
The singularity of g{t) for t —> 0 and the locality condition g{T) = 0 retrospectively justify the 
introduction of the cutoffs tscai and Tgcai for the temporal scaling regime. We could also have 
started from the spatial scaling relation m to obtain exactly the same functional form for g{t) 
with 

( 18 ) 

Thus the set of scaling relations d and m are compatible under the assumption of a constant 
weight-function /i = 1, i.e. there exists a solution for g(t) that satisfies (j1 411 and m simulta¬ 
neously. This sheds some light on the property of the weight-function h to select compatible 
temporal and spatial two-point correlations: scaling relations are among the simplest functional 
forms and allow h = 1, while for more advanced studies, such as deviations from scaling for 
At ^ [tscai, Tscai] and Ax ^ [heal, Lseal]. Other weight-functions h might be in order. For a brief 
account of this topic we refer the reader to M- 

To complete the specification of g{t), functional forms in the time intervals 0 < t < tscai and 
Tscai A i A T are needed in principle. For this analytical treatise, however, it is not necessary 
to specify the functional form of g{t) explicitly for these two time intervals, since we neglect 
the constants of proportionality Cx and ct in m and d in any case; examples addressing 
this issue can be found in Ref. d- The important point is that the validity of the scaling 
relations d) and d is independent of a specific choice of g{t) for t ^ [tscai, Tgcai]' Figure 1 
and Figure 2.a show that, for purely temporal separation, spatial scales in S larger than Ls^ai do 


obtain the 
(17) 


7 







not contribute to the ambit overlap Si fl S 2 for At > tscai, while spatial scales in S smaller than 
Iscai are completely part of the overlap for At < Tgcai and thus contribute only a term constant 
in At. 

Similar results hold for the purely spatial separation shown in Figure 2.b: regions of S smaller 
than I seal do not contribute to the overlap for Ax > I seal ■ The contributions from the large scales 
> Lscai result in a constant, as can easily be seen from 


F(Ax,0) = / ( 2 gr(t) — Ax) dt 

Jo 


C^scal 


2g{t)dt + 


.g(-l)(Ax/2) 


25 '(t)dt - 


r(2) 


K[2] - 2K[1] ’ 


(19) 


where denotes the inverse of g. Thus, a specific choice of g{t) for t < tgcai only involves 

the constant Cx and does not influence the scaling behaviour (d as such. The only restriction 
is 1 /( 0 , 0 ) < 00 (for finite expectations). 


3.3 Structure of higher-order correlations 

In the previous Section, we specified the model starting from scaling two-point correlations. It is 
now straightforward to derive scaling relations for all higher order correlations of purely spatial 
and temporal type. Section [3.41 shows how these scaling relations imply multifractality. 

First we note that, since h = 1, equation m translates to 


i(xi,ti;... ;xn,tn) = exp < / K 

( JRxR 




.2 = 1 


dx dt 


( 20 ) 


The argument of the cumulant function K in (PH) is piecewise constant where Yl^=i ^S{xi,u)i^j 0 
counts the number of field amplitudes e{xi,ti) that contribute to (x,t) via their ambit sets 
S{xi,ti). This function vanishes outside of UILi S{xi,ti). 

Focusing first on purely spatial two-point correlations of higher order, we get, using (PH), 
the analog to m 


(e(x, t)"'^e(x-b Ax, t)”^) = (e(x, (e(x-b Ax, t)"^) 

X exp {l/(Ax,0) (K[ni-b ^ 2 ] - K[ni] - K[n 2 ])} . (21) 

Translational invariance and d, d imply scaling relations for the higher-order two-point 
correlations 

(e(x, t)"'^e{x + Ax, oc (Ax)"^^"^!’"^)^ ^ Lscai], (22) 

where 

r(ni, na) = ^[2]'^-^^2K[l] • (23) 

(Again, the convexity of K implies K[ni -b n 2 ] — K[ni] — K[n 2 ] > 0.) The scaling range of (j22l) 
is identical to the scaling range of d and does not depend on the order (ni,n 2 ). 

An analogous procedure leads to scaling relations for the spatial higher-order three-point 
correlations illustrated in Figure 3. For ordered points xi < X 2 < X 3 with relative distances 
assumed to be within the spatial scaling range, \xi — Xj\ G [heal, Lgcai], (bJ = 1)2,3), we find 
that 


{e{x^,tr^e{x2,tre{xs,tD 

oc (x 2 -Xl)-"("l’” 2 )(x 3 -X 2 )■"^”^’”^)(x 3 -Xl)-«(”l’’" 2 ’" 3 ^ (24) 
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with a modified exponent ^ defined by 

^(ni, n 2 , na) = r(ni + n 2 , ns) - r(n 2 , ns) . (25) 

The reason for the different forms of the exponents r and ^ lies, of course, in the different ambit 
set overlaps: as shown in Figure 3, points xi and xs have only the one neighbour X 2 , while X 2 
has two. 

Equation (EU can be viewed as a generalised fusion rule in the sense of m- It is easily 
generalised to n-point correlations of arbitrary order because all overlapping ambit sets can be 
written as a combination of overlaps F(|xi — Xj|, 0 ) oc In |xj —Xi|, as long as |xj —Xj| € [IscahLgcai] 
for all point pairs. As shown by induction in m, the spatial n-point correlation for ordered 
points xi < X 2 < ■ ■ ■ < Xn and arbitrary order (mi,..., m^) satisfying — Xi € [Iscai, Lgcai] 
has the following structure: 

{eixi,tr^---e{xn,tr-) 

C n—1 \ n—1 n 

n (xi+i - n n , ( 26 ) 

i=l / j=2l=j+l 

where 

...,mi)= T{mi_j -b ... + - T(m/_j+i -b ... + m/_i, m/). (27) 

The modified scaling-exponents ^(mi,... ,mj) correspond to (j25|l for j = 3 and arise from the 
nested structure of the overlapping ambit sets. Physically, Eq. (OHl) implies that spatial n-point 
correlations factorise into contributions arising at the smallest scales Xj+i—Xj, at next-to-smallest 
scales Xj +2 — Xi, and so on up to the largest scale, x^ — xi. 

To complete the discussion of n-point correlations, we state the corresponding relation for 
temporal n-point correlations of arbitrary order 

{e{x,hr^---e{x,tnr'^) 


C n—1 \ n—1 n 

n iu+i - n n , (28) 

i=l / j=2l=j+l 

for ordered times ti < ... < tn and \ti — tj\ G [tscahTgcai], hJ = 1,... ,n. Finally it is to be 
noted that relations (PH) and (P|) only hold for purely spatial and purely temporal n-point 
correlations respectively. The general case of arbitrary n-point correlations (PH) does not allow 
a similar description in terms of scaling relations, since V{Ax,At) includes mixed terms in Ax 
and At. For a complete discussion of general space-time two-point correlations, we refer again 
to [Th) . 


3.4 Link to classical multifractality 


We complete the discussion of the multiscaling model with an investigation of the relation 
between multiscaling (defined as scaling of n-point correlations (1261) and (1281) 1 and classical 
multifractality (defined as scaling of coarse-grained moments). In the Appendix, we prove that 
multiscaling implies multifractality in the large scale limit. 

The term multifractality in the classical sense refers to n-th order moments of the field, 
coarse-grained at scale I centered on locations a, displaying scaling behaviour with some non¬ 
linear multifractal scaling exponent ^(n) > 0, 


Mn{cr,l) 


/1 j‘(T+l/2 

\ I Ja-l/2 





(29) 
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Note that this relation applies to stationary processes e{a) since the right hand side of (El is 
independent of the location a. Differentiating this relation twice with respect to I, it follows in 
second order, due to stationarity, that the two-point correlations 

{e{a + l)e{a)) oc (30) 


scale with the same scaling exponent ^(2) as M 2 . The inverse need not be true, for the scaling 
relation El becomes singular for 1^0, though at small scales deviations from El have to 
occur which in turn may destroy the relation El HlEl. However, El indicates a strong 
connection between scaling of n-point correlations and multifractal scaling of order n. 

The multiscaling model implies scaling relations for n-point correlations, with deviations 
from pure scaling for scales smaller than Igcai for spatial correlations and scales smaller than 
tscai for temporal ones. Thus the question arises whether multifractal exponents n{n) are to be 
expected (see also m ). To answer this question, we assume the one-point moments {e{x,t)^) 
to be finite (i.e. we restrict to Levy bases with K[n] < 00 ). 

In the Appendix it is shown that the integral moments of temporal type 



asymptotically exhibit scaling behaviour for tgcai ^ L Moreover this is also true for the integral 
moments of spatial type 



for heal ^ I with the same multifractal scaling exponents 


/i(n) 


t(2) 

K[2] - 2K[1] 


(K[n] - nK[l]). 


The crucial assumption that enters the proof of m and El is 


(33) 


.o^ K[n]-K[n-l]-K[l] 
^ ’ K[2] - 2K[1] 


/i(n) — fi{n — 1) < 1. 


(34) 


This assumption ensures that large scale correlations dominate the moments of the coarse grained 
field. It is to note that El is a sufficient condition for multifractality in the large-scale limit. 
The statistics of the energy dissipation in fully developed turbulence is an important example 
of an observable where condition El holds; see Ref. mi 

The identity of spatial and temporal multifractal scaling exponents n{n) is clearly a result of 
the identical scaling behaviour of purely spatial and temporal n-point correlations. The scaling 
of spatial and temporal integral moments is independent of the choice of the boundary function 
g(t) for t ^ [tscahTgcai] as long as 14(0,0) < 00 . Under these mild restrictions, we are able to 
model a wide range of scaling exponents ^(n) by choosing a proper cumulant function K via the 
Levy basis Z that fulfils the sufficient condition EH)- Examples are ^(n) oc n“ — n for a stable 
basis with index of stability 0 <a< 2 ,a 7^1 and 


g,{n) oc (1 — n)\J— ( 5 "^ + n\J— {[5 + 1)^ — \J — {(5 + n)^, (35) 

with \(5 + n\ < a, for a normal-inverse-Gaussian distribution NIG(q;,/?, d, i/) EHHSim]- De¬ 
pending on the parameters that characterize the distributions, there exists a critical order ric 
where EU does not hold any more. The NIG(a, /3, S, u) distribution is an example of a Levy 
basis where multifractality EHI) is defined only up to a finite order n, since K[n] < 00 only for 
\P + n\ < a; for larger n, the moments (e”) and M„ do not exist. 
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4 Conclusion 


We have presented a general framework for modeling of spatio-temporal processes that allows, 
even in its generality, an analytical treatment of general spatio-temporal n-point correlations. 
This framework consists of a homogeneous Levy basis, the concept of an ambit set as an asso¬ 
ciated influence domain and a weight-function h. These three degrees of freedom can be chosen 
arbitrarily and independently, thus encompassing a wide range of applications. In this respect, 
we mentioned briefly related work [niiiiiis] and showed them to be special cases of this frame¬ 
work. In a specific illustration, we have shown that a stationary and translationally invariant 
version of the general model can be used to construct a multiscaling and multifractal causal 
spatio-temporal process starting from scaling relations of two-point correlations. 

Many applications immediately come to mind. The great flexibility and tractability of the 
framework’s mathematics might well find its way into modeling of rainfields, cloud distributions 
and various growth models, to name just a few examples of spatio-temporal processes. Another 
field of application for the special case of the multiscaling model is the description and modeling 
of the statistics of the energy-dissipation in fully developed turbulence as a prototype of a multi¬ 
fractal and multiscaling field. A first step in this direction was undertaken in m, where scaling 
two- and three-point correlations m and (EHl) were shown to be in excellent correspondence 
with data extracted from a turbulent shear flow experiment. 


Acknowledgements 

This work was supported in part by the EU training research network DynStoch - Statisti¬ 
cal methods for Dynamical Stochastic Models. O.E.B.-N. and J.S. acknowledge support from 
MaPhySto - A Network in Mathematical Physics and Stochastics, funded by the Danish Na¬ 
tional Research Foundation and support from the Alexander von Humboldt Foundation with 
a Feodor-Lynen-Fellowship. H.C.E. acknowledges support from the South African National 
Research Foundation. 


Appendix: Scaling relations for integral moments 


This appendix proves the classical multifractal property (PH) for the multiscaling model in the 
limit tscai 'C t < Tgcai and Igcai ^ ^ < Lg^ai under the assumption that 


-K[2]-2K[1]- < 


(36) 


The proof is carried out in detail only for the spatial case; the temporal counter part of the 
above statement is straightforward. 

With the abbreviation dnih, ■ ■ ■ ,ln) = ■ ■ ■ e{ln,t)) , 0 < ^2 < • • • for the 

spatial correlation function of order n and using the translational invariance of the correlation 
structure, the spatial integral moments of order n (El are given by 


1) = n\ f din t dln-l ■■■ t dh (/ - In) dn (h, ■ ■ ■ , L)- (37) 

JO Jo Jo 

To calculate the involved overlaps of the influence domains, it must be distinguished whether the 
spatial distances are smaller or larger than I seal- In the limit I seal ^ I, the dominant contribution 
is 


mW(x,o^mW(0 


(38) 


II 



r^n ^scal r^3 ^scal 

— 7z! ^ I dl-n—X * * * / ^^2 (^ ^n) ^n(^2? • • • ; ^n)* 

l)^scai d (^n 2)lscal ^ Ig^al 

The proof of the multifractality of Mn'^ {x, 1) is carried out in two steps. The first part shows that 
Mn\l) oc in the large scale limit. In the second step, we show that the approximation 

Ain {x,l) ~ Mn {1) holds for I ^ heal- We also provide a rough estimate for the relative error 

The correlation function dn can be rewritten with the help of the generalised fusion rules 
(12(11) as 


n k—1 

k=2j=l 

where 

times 

and ^2 = ^(l, 1). In the next step, we define 

Pn {^ 1 1 seal) 


(39) 

(40) 

(41) 


= r 


fl 

pin Iscal 

r 

1 dlfi 

1 dln—l ■ ■ 

• / 

l)^sca/ ‘ 


dls 


^3 ^scal 


n k—1 

dh (I - ^n) n n 

k=2j=l 


-^j+l 


Note that Mn{l) oc Fn{l,lscai) with a constant of proportionality that is independent of 1. 
With the abbreviation 


it follows from (US and (f?7|) that 


where 


k—1 

h{k) = -j2^,+i, 

(42) 

i=i 


n 

h{k) = -Kn) 

(43) 

k=2 


> ''^1 k[2] - 2K[1] ■ 

(44) 


Thus we get, using condition (jsni) 


. ,KM -K[n- 11 - Km 

h{n) = fi{n - 1 ) - fJ-{n) = -t( 2 )— —rvTTri - ^ > — 1 - 


K[2] - 2K[1] 


(45) 


It follows immediately that 


Mnil)P^^^ OcFnil,lscal/l)- (46) 

Fn0-,lscai/^) is positive and increasing with increasing 1. It is easy to show that Fn{l,lscai/l) is 

jh{ 

‘‘k 

1 


bounded. From dsi and dH it follows that dn < m=2 and therefore 


/ I p1 p1 ^ ^ 

din / df„_i ... d/2n^ n 

-scal/^ '^^sca l/^ '^^scal/^ fc = 2 k = l 


_i 1 + h{k) 


(47) 


The last step in (|17|) requires (|HH|) to hold. Since Fn{l,lscai/l) is increasing with I and bounded, 
there exists a constant c with 

lim M„(/)/^(’"i = c < oo (48) 

l—*^oo 
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and therefore 


Mn{l) OC 


(49) 


in the large scale limit I » I seal- 

To complete the calculations, we provide a rough estimate of the relative error between the 
exact relation and its approximation PS|). By going from to we neglect all n-point 
correlations with one or more distances \li — Ij\ < heal- These are (^) integrals of the form 


nir^ f din--- 

iP' l)^sca/ 


r^i-\-2 ^scal 


dk 


i+l 


H+l 




dh 




'( 2 - 2)4 


d^ 2—1 • • • 


f's 4 


^scal 


(50) 


dl 2 {l ^n)^n(^ 2 ) • • • ) ^71)5 


where one distance (chosen to be k+i — k in (1^171) 1 is smaller I seal and ( 2 ) integrals where two 
distances are simultaneously smaller I seal etc., and one integral where all distances are smaller 
than Igcai- Each of these integrals have an upper bound ..., 0) (assumed to be 

finite), where k denotes the number of distances that are smaller I seal- Thus we have 


( x , 0 - ( x , 0 < nll-^J2u] tar~''dn ( 0 ,... , 0 ) 


k=l 


= nU-^dn{0,...,0){{lscal + ir-n. 


The relative error 






< const. Xl^^^'>-^{{lseal + ir-n 


(51) 


(52) 


tends to zero for I 00 and n > /i(n) (which is always true, for l"'Mn\x,l) is monotonically 
increasing for positive and finite n-point correlations). The results (1491) and isi are independent 
of the choice of the small scale statistics as long as they are finite. 
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e{x, t] 



X + git' — t + T) 


Tscal 


Figure 1; Illustration of the spatio-temporal ambit set S{x,t) associated with the field amplitude 
e{x, t) and bounded by a monotone function g{t' — t + T). 


(x, t + At) 


e{x + Ax, t + At) 


e{x,t) 


e{x,t), ,e{x + Ax,t) e{x,t) 


m 


Figure 2: Spatio-temporal overlaps (shaded areas) of the ambit sets separated by a (a) temporal 
distance At, (b) spatial distance Ax and (c) spatio-temporal distance {Ax, At). 















